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AUTHOR* S  SUMMARY 

The  collision  of  a  spherical  liquid  drop  with  a  plane  rigid  surface  at  a 
right  angle  is  considered  and  an  equation  derived  for  determining  the  pressure 
in  the  zone  of  contact.  The  results  of  experimental  determination  of  the 
pressure  are  given  which  allow  the  theoretical  equation  obtained  to  be  used 
for  pressure  determination  in  the  solution  of  practical  problems. 
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In  estimating  the  wear  of  surfaces  of  aircraft  due  to  rain  and  wear  of 
turbine  blades  by  moist  steam  we  use  an  equation  to  determine  the  maximum 
pressure  arising  at  the  point  of  impact  of  a  liquid  drop  against  a  plane  solid 
surface.  Such  an  equation,  proposed  in  Ref.l,  is  also  given  in  Refs. 2  and  3, 
i.e. 


P  -  fpVc^ij.  (I) 

where  a  is  a  coefficient  characterising  the  rigidity  of  the  drop,  p  is  the 
normal  density  of  the  liquid,  V  is  the  velocity  of  water  particles  behind  the 
shock  wave  taken  as  equal  to  the  normal  collision  velocity,  c^j^  is  the  velocity 
of  propagation  of  the  front  of  the  shock  wave  in  the  liquid  which  increases 
with  the  rise  in  pressure  in  accordance  with  the  data^.  However,  experimental 
findings^* ^  show  that  in  the  range  of  collision  velocities  150-600  m/ti  the 
pressure  at  the  point  of  impact  corresponds  to  the  equation 

P  -  pVc^j^.  (2) 

Thus,  experimentally  determined  pressures  at  high  velocities  correspond 
to  a  value  of  2  for  a  in  equation  (1)  which  is  contrary  to  the  physical  sense 
of  the  coefficient  a  which  cannot  be  greater  than  unity*,  and,  on  the  other 
hand,  equation  (2)  does  not  agree  with  experimertal  findings  for  collision 
velocities  below  150  m/s. 

Let  us  consider  the  collision  of  a  fixed  spherical  drop  with  a  plane 
surface  with  an  approach  velocity  V.  In  Fig. la  the  position  of  the  drop  and 
the  plane  is  shown  at  the  moment  of  impact  (t^  <*  0)  when  the  plane  is 
projected  on  the  x-axis  of  the  system  of  co-ordinates  xoz  and  the  lower  part 
of  the  drop  is  represented  by  the  arc  BOC. 

The  line  1-1  is  the  projection  of  the  same  plane  at  the  moment  tj. 
Experiments  show  that  the  drop  at  this  moment  consists  of  a  spherical  undeformed 
part  and  a  cylindrical  part  BDEC  separated  from  it  by  the  X-plane  and  formed  as 
a  result  of  two  simultaneously  occurring  deformations  -  the  elastic  deformation 
of  part  of  the  drop  from  position  BOC  to  position  BAG,  and  the  deformation 
resulting  from  the  expression  of  the  liquid  from  position  BAC  to  position  BDEC. 
It  is  considered  that  BAC  is  a  parabola,  and  insofar  as  the  oppositely  hatched 
parts  correspond  to  equal  volumes  of  liquid  BD  >  AF  °  6,  as  shown  in  Fig.  la. 
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Following  the  instant  -  0  begins  the  penetration  of  the  plane  into 
the  drop,  and  this  is  characterised  at  the  instant  by  the  segmen'  AF  -  5 
and  by  translation  of  the  plane  for  a  distance  £'  >  £.  For  the  time  interval 
At  ■  tj  ~  tQ  the  relations  obtain, 

6*6  6 

At  -  —  "  Tj"  •  ^1  “  F  » 


where  Vj  is  the  velocity  of  motion  of  the  plane  with  respect  to  the  liquid. 
For  the  absolute  velocity  of  the  liquid  v  we  obtain  the  expression 


where 


V 


V  -  aV  , 


(A) 


a 


(5) 


For  velocity  Vj 


Vj  -  V  -  V  -  (I  -  a)V 


(6) 


The  velocity  v  of  motion  of  the  X-plane  in  the  liquid  is  twice  as  high  as  the 
velocity  Vj  since  in  time  At  it  is  translated  within  the  liquid  for  the 
distance  26 

V  -  2(1  -  a)V  .  (7) 


The  mass  of  liquid  drawn  into  movement  with  velocity  oV  is  given  by 


m  ■  fflj  +  m2 


(8) 


where  mj  ■  irR  6p  is  the  mass  of  liquid  drawn  into  motion  below  the  A-plane 
(in  turn,  R  ■  FD  ■  FE),  and  mj  ■  is  the  mass  of  liquid  above  the 


A-plane  (Fig.  lb)  drawn  into  collision  by  the  shock  wave  of  velocity  c  .  passing 

sn 

through  the  drop.  Assume  that  the  volume  bounded  by  the  parabolic  surface^ the 
projection  of  which  on  the  plane  xoz  is  shown  by  the  dotted  linek  is  equal  to 
the  volume  of  a  cylinder  of  half  the  height  and  with  the  original  base.  To 
determine  the  mass  m2  two  states  at  times  tj  and  tj  -  A't  are  taken. 
Correspondingly,  zA  and  z  are  the  co-ordinates  of  the  A-plane  at  times 
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tj  and  tj  -  A*t;  A*z  is  the  distance  travelled  by  the  X-plane  in  time  A*t; 
z*  is  the  unknown  ordinate  of  an  arbitrary  point  of  the  surface  bounding  the 
mass  of  the  liquid  drawn  into  motion  above  the  X-plane;  and  z^^  is  the 
maximum  unknown  ordinate  of  this  surface  at  time  tj.  Then  for  the  time 
interval  A't  the  following  relations  obtain 

. , .  A'z  +  z'  A'z 
A't  ■  — — -  ■  — —  . 


whence 


(9) 


Clearly,  A’z  -  Zj^  -  z,  where  Zj^  -  6'  +  6  (Fig.l).  Taking  equation  (5)  into 
account,  after  simple  transformations 

1  -  - 

We  should  note  that  the  expression  for  in  Ref.  I  does  not  contain  the 
1  -  2. 

term  - -  and  it  is  assumed  without  explanations  that  2,-26.  This  excludes 

1-0  ^ 
consideration  of  the  elastic  deformation  of  the  drop  when  the  ordinate  z^  is 

being  determined.  The  consequence  of  this  assumption  is  increasing  error  in  the 

determination  of  z^  as  the  part  played  by  elastic  deformation  increases  with 

increasing  collision  velocity.  Furthermore,  as  the  point  M  at  time  tj  -  A't 

is  part  of  a  sphere,  its  ordinate  z  is  determined  from  the  equation  of  a 
2  2  2  2 

sphere  r  -  x  +  y  +  (r  -  z)  ,  and  hence  when  y  «  0 

1  -  1 

n  9  I  *2  2  2  i 

z  ■  r  -  (r^  -  and  A*  z  ■  Zj^  -  z  ■  26  y  -  -  r+(r  -x)*. 


Using  equation  (9)  we  can  write 


and  then 


z 

max 


-  26 


1  -  a 


-  26 


(10) 
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As  a  result  of  substituting  the  value  of  z*  in  the  expression  for  m. 

max  2 

1  -  - 

m2  -  A— I  ,  (II) 

and  after  substitution  of  the  values  for  mj  and  m2  in  the  equation  for  m 

1  -  2. 

2  2*2 

m  »  fflj  +  m2  ■  irpSR  +  irpfiR  $  ^  .  (12) 


We  now  make  use  of  the  momentum  equation 


mAv 

p  ■ 

At 


» 


(13) 


where  Av  ■  v(tj)  -  v(tQ)  ■  oV,  since  v(tQ)  ■  0. 

For  the  pressure  arising  in  the  collision  of  a  drop  with  a  plane  we  get 


P  mV 

2  2 
ttR  AtnR^ 


(14) 


Taking  into  account  the  fact  that  At  -  _  p)y  »  after  substitution  of  the 

value  of  m  from  equation  (12)  and  appropriate  transformations 


o  ■  2  V  ] 


(15) 


Thus  in  contradistinction  to  equation  (I),  equation  (15)  contains  the 
V  \ 

-  a  - — I  and,  for  example,  even  when  V/c  .  ■  1,  when  a  ■  0.8, 
%h/ 

the  ratio  P(i5)/P(])  "  2.2. 

To  make  calculations  using  equation  (15)  it  is  necessary  to  find  values 

for  a  or  k  over  a  wide  range  of  velocities.  Treatment  of  experiments^ 

8 

taking  account  of  an  analytical  study  allows  us  to  conclude  that  the  value  of 

V 

k  reaches  unity  when  —  ^  1  corresponding  to  the  velocity  of  sound  in  water 

c  ■  1500  m/s.  The  experimental  data  in  Ref. 6  confirm  this  conclusion.  One 

•  V  “3 

experimental  value,  a  ■  0.4  when  -  7.3x10  is  given  in  Ref.l.  The 

c 

characteristic  calculation  point  corresponding  to  the  intersection  of  the 
0  and  k  curves  in  equation  (15),  i.e.  where  a  ■  k,  corresponds  to  a  ■  0.67. 
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We  carried  out  experiments  to  elucidate  further  the  unknown  relationships. 
At  collision  velocities  up  to  V  >  6  m/s  drops  were  hurled  at  a  plane  surface 
of  barium  titanate  piezoceramic  crystal  of  diameter  20  mm  and  thickness  1.5  mm 
mounted  in  a  special  'anvil'  (Fig. 2).  The  dropper  and  liquid  were  located 
at  heights  ^1234*  1.15,  0.3  and  0.25  metres.  The  diameter  of  the 

drops  was  3  mm.  Their  impact  at  a  constant  point  of  the  crystal  was  assured 
by  a  cone  mounted  above  the  crystal  with  an  orifice  of  diameter  d  >  4  mm  at  the 
apex.  The  signal  from  a  transducer  was  transmitted  to  a  pulse  oscillograph 
mark  S-1-19,  and  a  'Zenith'  camera  with  in  attachment  was  used  to  take  films 
from  the  screen.  The  piezoelectric  properties  of  the  barium  titanate  ceramic 
subjected  to  high-voltage  polarisation  at  low  mechanical  loads  (up  to 
1800  kg/cm  at  15  C)  are  linear  .  The  linearity  of  the  characteristics  of  the 
crystal -oscillograph  system  was  confirmed  by  special  experiments.  Due  to  the 
fact  that  drops  of  constant  diameter  were  of  different  rigidity  depending  on 
the  velocity  of  collision  the  direct  proportionality  of  the  height  of  the  blips 
on  the  screen  to  the  impact  velocity  was  impaired,  and  this  made  it  possible, 
using  equation  (15)  to  estimate  the  change  of  a  in  comparison  with  its  known 
value  for  the  maximum  collision  velocity  and  to  obtain  a  number  of  values  of  a. 

Another  method  of  determining  a  consisted  in  hurling  drops  at  the 

polished  surface  of  a  copper  disc.  The  drop  in  flight  passed  through  the 

beam  of  a  photoelectric  cell  which  after  a  time  lag  to  allow  the  drop  to  reach 

3 

the  target  switched  on  an  SKS-1  camera  operating  at  4  x  10  frames  per  second. 

When  the  films  were  processed  it  was  possible  to  determine  the  velocity  of 

translation  of  the  plane  through  the  drop  when  the  collision  velocity  was 

5.7  m/s,  and  from  equation  (6)  to  get  the  value  a  ■  0.3.  Finally,  an  SSKS-1 

4 

camera  operating  at  6  x  lO  frames  per  second  was  used  to  film  the  collision 
of  a  drop  with  a  target  rotated  by  an  electric  motor  (Fig. 2)  at  V  -  30  m/s. 

As  a  result  of  processing  five  films  it  was  possible  to  determine  the  velocity 
of  translation  of  the  plane  through  the  drop,  a  being  found  to  be  0.6. 
Generalisation  of  available  experimental  data  and  of  our  own  experiments  led 
to  our  obtaining  the  dependences  of  a  and  k  on  V  (Fig. 3)  to  be  used  in 
determining  pressures  from  equation  (15),  which  in  contradistinction  to 
equations  (1)  and  (2)  is  applicable  for  normal  collision  velocities  from 
0  to  600  m/s.  At  higher  velocities,  the  expenditure  of  part  of  the  collision 
energy  on  deformation  of  the  surface  and  the  drop  leads  to  a  certain  reduction 
of  pressure  in  comparison  with  calculated  (theoretical)  values.  The  pressures 
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calculated  from  equation  (15)  for  a  number  of  practical  cases  allow  of  a 
satisfactory  explanation  for  the  observed  intensive  erosion  damage  of  surfaces 
of  aircraft  components  coming  into  collision  with  raindrops  and  of  turbine 
blades  eroded  by  moist  steam. 
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Fig.] 
Fig.  2 


Fig.  3 


Diagram  for  determining  the  mass  of  a  drop  drawn  into  collision: 
a  -  below  the  X-plane;  b  -  above  the  X-plane. 

Apparatus  for  the  study  of  the  collision  of  a  drop  with  a  solid 
surface:  1  -  frame;  2  -  cone;  3  -  drop  catcher;  4  -  glass  tube; 

5  -  dropper;  6  -  electric  motor;  7  -  moving  target;  8  -  cooled  target; 
9  piezocrystal  target;  a  -  crystal;  b  -  body;  c,  d  -  insulators; 
e  -  contact  diac;  f  “  spring;  g  -  electric  terminals;  10  -  camera; 

II  -  light  source;  12  -  drop.  Below  right:  diagram  of  method  of 
filming  the  process  of  collision  of  the  drop  with  the  target  and 
electronic  commutating  circuit:  I  -  light  source  (27  volt  lamp); 

II  -  photocell  type  STs8-3;  III  -  cathode  follower;  IV  -  amplifier; 

V  -  electronic  relay;  VI  -  rectifier. 

Experimental  dependences  of  the  coefficients  k  and  a  on  V,  plotted 
on  the  basis  of  generalised  experimental  datai  I  -  k,  2  -  a. 
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